Crystal-field and Kondo scale investigation of CeMIns (M=Co, Ir and Rh): 
a combined x-ray absorption and inelastic neutron study 
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Linear polarized soft-x ray absorption (XAS) and inelastic neutron scattering (INS) experiments 
have been performed on CeMIns with M = Rh, Ir, and Co to determine the crystal-field scheme and 
characteristic Kondo temperatures T* for the hybridization between 4/ and conduction electrons. 
The ground state wave functions are determined from the polarization dependent soft-XAS data 
at the cerium M4,5 edge and the crystal-field splittings from INS. The characteristic temperature 
T* has been determined from the line widths of the neutron scattering data. We find that the 
quasielastic line widths of the superconducting compounds CeCoIns and Celrlns are comparable 
with the low energy crystal-field splitting. 

PACS numbers: 71.27.+a, 75.10.Dg, 78. 70. Dm, 78.70.Nx 
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I. INTRODUCTION 

The ternary rare earth family CeMIns (M = 
Co, Ir, and Rh) are heavy fermion, unconventional 
superconductors^^ CeCoIns and Celrlns become super- 
conducting at ambient pressure at Tc — 2.2^ and Tc = 0.4 
K"^, whereas the incommensurate heavy fermion antifer- 
romagnet CeRhlns (Tn = 3.8) exhibits superconductiv- 
ity under pressure {Pc = 1.6 GPa, = 2.1K)^ AU mem- 
bers of this family crystallize in the tetragonal HoCoGas 
structure (space group P4/mmm), which is derived from 
cubic Celn3 intercalated with MIn2 layers along the 
tetragonal c axis. The cubic compound Celns orders an- 
tiferromagnetically at Tjv = 10 K and has a hybridization 
temperature T* of about 10 K.^ This temperature, which 
gives the energy scale of the hybridization between the 
local 4/ moments and surrounding conduction electrons, 
is highly pressure dependent in CeMIns compounds. It 
has been argued that a sufficiently strong hybridization 
or Kondo interaction with respect to the RKKY ex- 
change interaction suppresses antiferromagnetic order to 
the benefit of superconductivity.^'^ General scaling be- 
haviors of the characteristic energy scales in Kondo lat- 
tice materials are a matter of intense debate^'"— and 
T* has traditionally been interpreted as the temperature 
below which the coherence of the Kondo lattice sets in. 
More recently T* has been suggested to be the temper- 
ature scale denoting the development of a dense " Kondo 
liquid" within a "two- fluid" model also comprised of a 
"Kondo gas" phase of uncorrelated magnetic moments 
at high temperatures. This characteristic temperature 



scale is determined with many techniques such as ther- 
modynamic, transport, knight shift measurements, and 
of course, quasielastic neutron scattering^^"^ Another 
interesting scaling has been shown by Bauer et ali^: 
the superconducting transition temperatures Tc of the 
CeMIns and also of the PuMGas family vary linearly 
with the c/a ratio of the tetragonal lattice constants, 
pointing towards the importance of the anisotropic elec- 
tronic structure for the superconducting state. This 
brings into focus the importance of the spatial distribu- 
tion of the crystal-field split Hund's rule ground state, 
which is highly anisotropic for materials containing rare 
earth. 

The Hund's rule ground state of Ce'^+ with J = 5/2 
splits under the influence of a tetragonal crystal-field 
(point group D4h) into three Kramer's doublets, which 
can be represented in the basis of \Jz)- The eigenfunc- 
tions of the three Kramer's doublets can be written as 



|2) 
|1) 
|0) 



= Fe = I ± 1/2) 
= T\= /3|±5/2)-a|T3/2) 
«|±5/2)-H/3|t3/2) 



(1) 



F2 
J- 7 



with + ji^ = 1. The anisotropy of certain wave 
functions may give rise to 4/ conduction electron hy- 
bridizations which are more advantageous than others for 
forming a superconducting ground state. The impor- 
tance of momentum dependent (q-dependent) hybridiza- 
tion in these and some semiconducting Kondo materials 
has been investigated by several groups.—"— 

Various groups attempted to determine the crystal- 
field scheme of these compounds, but there are significant 
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discrepancies depending on the applied methods, which 
include bulk measurements based on transport, ther- 
modynamic and NMR experiments.— i^"— Christianson 
et al. performed extensive inelastic neutron scattering 
(INS) studies,—!^ but phonon contributions in the en- 
ergy range of the magnetic scattering, broadened crystal- 
field excitations due to hybridization effects and the enor- 
mous absorption of the sample's constituents make the 
determination of reliable magnetic intensities rather chal- 
lenging. Since the latter give the wave functions via the 
transition matrix elements the resulting wave functions 
should be taken with care (see appendix) while the tran- 
sition energies are fairly sound (see section III B). 

We have shown for the case of the heavy fermion mate- 
rials CePd2Si2 and CePtaSi that polarization dependent 
soft x-ray absorption (XAS) at the Ce M4_5 edges can be 
complementary to neutron scattering in determining the 
ground state wave function. '^^'^^ XAS is highly sensitive 
to the initial state and via its polarization dependence 
(linear dichroism LD) direct information about the \ Jz) 
admixtures of the ground state wave function can be ob- 
tained. Sensitivity to higher lying crystal-field states is 
achieved by thermally populating those statesi '^^'^^ 

We present a combined inelastic neutron scattering and 
soft x-ray absorption study on the CeMIus M = Co, Ir, 
and Rh compounds. Combining both techniques has the 
advantage of determining transition energies and mixing 
parameters independently, each with the most suitable 
technique. The line positions in the magnetic contribu- 
tions of the INS data yield the crystal-field transition 
energies within meV resolution, whereas the LD in XAS, 
when performed at temperatures where only the ground 
state is populated, yields the ground state wave function, 
i.e. in case of a mixed ground state the mixing factor 
(see eq. (1)) within Aa^ — ±0.0025. In the limit of small 
crystal-field splittings the dichroic signal of the ground 
state is independent of crystal-field energies. Once the 
ground state has been determined, the order of states 
can be determined from the temperature dependence of 
the LD since at finite temperatures it refiects the su- 
perposition of polarizations from each populated state, 
i.e. here the crystal-field energies, as determined from 
the neutron scattering experiment, enter via the ther- 
mal population. In addition, high-resolution INS data 
are presented to determine the characteristic tempera- 
ture T* for the 4f conduction electron hybridization via 
the quasielastic line width. The latter has been applied 
successfully by several authorSi ^^'^'^" — 



II. EXPERIMENTAL AND DATA 
CORRECTION 

The high quality single crystals of CeMIus for the x-ray 
experiments were grown with the fiux-growth method^. 
The powder samples for the present INS experiments 
were the same samples as used by Christianson et al^^^. 



A. XAS 

The XAS spectra were recorded during various beam 
times at the two different synchrotron light sources 
BESSY II and NSRRC. We recorded ah spectra with the 
total electron yield method (TEY) and under UHV, i.e. 
in a chamber with a pressure in the 10"^° mbar range. 
Clean sample surfaces were obtained by cleaving the sam- 
ples in situ. At BESSY II we used the UE46 PGM-1 un- 
dulator beam line. The total electron yield (TEY) signal 
was normalized to the incoming photon fiux Iq as mea- 
sured at the refocusing mirror. The energy resolution at 
the cerium M4_5 edges {hi' w 875 — 910 eV) was set to 
0.15 eV. The undulator combined with a normal incident 
measurement geometry allow for a change of polarization 
without changing the probed spot on the sample surface. 
The two polarizations were E J-c and -E || c, c being the 
long tetragonal axis. At the NSRRC we performed the 
experiment at the Dragon dipole beam line. The energy 
resolution at the cerium M4_5 edges was set to 0.4 eV. 
The crystals were mounted with the c-axis perpendicu- 
lar to the Poynting vector of the light. By rotating the 
sample around this Poynting vector, the polarization of 
the electric field can be varied continuously from iJ _L c 
to E\\c. For all measurements the sample was rotated 4 
times by 90°, so that for each orientation E J-c and E\\c 
two equivalent positions were measured. Thus for both 
experimental end stations a reliable comparison of the 
spectral line shapes is guaranteed. We measured several 
crystals and/or recleaved in order to assure the repro- 
ducibility of the spectra (see table 1). 

To calculate the XAS spectra we performed ionic full 
multiplet calculations using the XTLS 8.3 program^-. 
All atomic parameters are given by Hartree-Fock val- 
ues, with a reduction of about 40% for the 4/ — 4/ 
Coulomb interactions and about 20% for the 3d — 4/ 
interactions to reproduce best the experimental isotropic 
spectra, /isotropic — 2/_l + /|| . These values compare well 
with our findings for other heavy fermion compounds^i^ 
and account for the configuration interaction effects not 
included in the Hartree-Fock scheme. Once the atomic 
parameters are fine tuned to the isotropic spectra, the 
polarized XAS data can be described by the incoherent 
sums of the respective polarization dependent spectra of 
the pure jj^) states^ as long as the crystal-field split- 
ting EcF is small with respect to the spin orbit splitting 
Eso- The latter requirement, which is fulfilled here (Eso 
Ki 280 meV and Ecf < 30 meV), assures that interfer- 
ence terms resulting from intermixing of the J — 5/2 and 
J = 7/2 multiplet can be neglected. 

B. Inelastic neutron scattering 

We have measured the neutron scattering function 
S{Q,uj) of CeRhlus, Celrlus, and CeCoIus with the in- 
elastic time-of-fiight spectrometer HET at the neutron 
spallation source ISIS with incoming energies of 20 and 
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FIG. 1: (color online) Low temperature (T = 20 K) linear 
polarized XAS spectra of CeRhlns, Celrlns, and CeCoIns at 
the Ce^^ M4,5 edges. The sohd lines are the measured data, 
the dotted ones the simulations as described in the text. The 
orbitals represent the spatial distribution of the 4/ wave func- 
tions according to the respective ground state admixtures, 
a| ± 5/2) + /3| T 3/2) 



transfers at elastic position for the low and high angle 
groupings are Q « l.SA and « 10. OA for the 60 meV and 
Q K, 1.0 A and ~ 5. 5 A for the 20 meV data. Some data 
were taken at the cold time-of-flight spectrometer IN6 at 
ILL with an incoming energy of 3.1 meV and an energy 
resolution of 70 /xeV at elastic position. All detectors 
from 10° to 115° have been grouped together. Because 
of the small incident energy the momentum transfer Q 
is <2 A in the energy window of interest. A flat sample 
geometry was used for all neutron experiments and - be- 
cause of the enormous absorption of In, Rh, and Ir - well 
defined, but thin samples were crucial in order to guar- 
antee transmissions of at least 30%. This reduced the 
sample amount for the IN6 experiment to about 5 g. All 
data have been normalized to monitor count rate and 
vanadium and have been corrected for absorption and 
self-shielding. The description of the phonon correction 
has been moved to the appendix. 



III. RESULTS 

For Ce'^"*" in T>Ah point symmetry the crystal-field 
Hamiltonian Hcf ^ B^O^ + 3^01 + BjOj describes the 
crystal field potential when the three Stevens parameters 
B2, B2, and Bf are determined. The BJP" parameters 
are determined via the crystal-field transition energies 
within the Hund's rule ground state and the mixing pa- 
rameter a. In section A the mixing parameter will be 
determined from the LD in the XAS data at low tem- 
perature, in section B the crystal-field energies are ob- 
tained from thermal neutron scattering data, and C the 
order of crystal-field states will be confirmed from the 
temperature of the LD effect in the XAS spectra. In sec- 
tion D we finally present the quasielastic results from the 
cold neutron data in order to determine the hybridization 
temperatures T*. 



sample synchrotron facility sequence of temperatures [K] 



CeCoIns #1 NSRRC 


78,12, 295 


CeCoIns #2 NSRRC 


10 


CeCoIns #3 BESSY 


18, 50, 80, 130, 17, 180, 280, 




17, recleave, 19, 50, 80, 280 


Celrlns #1 NSRRC 


18, 50, 80, 150, 300, 




recleave, 300 


CeRhlns #2 NSRRC 


20, 80, 295 


CeRhlns #2 BESSY 


18, 50, 80, 150, 300 



TABLE I: Experimental details about the XAS measurements 



60 meV and energy resolutions of 0.6 and 1.8 meV in the 
2.5 m detector banks. The low angle banks cover 26= 9° 
to 29° and the high angle banks 130° to 140°. All low and 
high angle detectors are grouped together respectively in 
order to gain statistics. We therefore refer to S{29,uj) 
from now on. The corresponding averaged momentum 



A. Low temperature polarized soft XAS: the 
ground state wave function 

Figure [T] shows the low-temperature 20 K linear polar- 
ized XAS data of CeRhlns, Celrlns, CeCoIns at the Ce^+ 
M4,5 edge full lines), i.e. at a temperature sufficiently 
low so that the ground state is populated. The latter 
has been verified with the knowledge of the crystal-field 
energies (see table 2). For CeCoIns the linear dichro- 
ism (LD) is smallest, but has an opposite sign with re- 
spect to Celrlns and CeRhlns. CeRhlns has the largest 
LD. Figure [1] shows further the simulations based on a 
full multiplet treatment as described in the experimental 
section (dotted lines) . For CeCoIns we find from our full 
multiplet calculation a = 0.36. This value is just below 

the zero polarization for a = « 0.41 where the LD 

vanishes when it changes sign. For Celrlns and CeRhlns 
we obtain a = 0.50 and a = 0.62, respectively. The re- 
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FIG. 2: (color online) Inelastic neutron data at 5 K with 20 meV a)-e) and 60 meV f)-h) incident energy. The blue triangles 
represent the total scattering, the black lines the phonon intensities as obtained from a high to low angle scaling. The red 
circles show the magnetic scattering obtained after phonon correction (see appendix) and the red lines are Lorentzian fits. 



[meV] CeRhlus Celrlus CeCoIns instrument CeAu2Si2 

r,„/2@ 8K 1.2 ± 0.2 IN6 3.1 meV OlT 

r,„/2Q 75K 1.7 lb 0.2 2.7 ± 0.5 3.9 ± 0.5 IN6 3.1 meV 0.43 

E[ 7.0 ± 0.5 5.2 ± 1.0 6.8 ± 1.0 HET 20meV 16^ 

E2 24.7 ± 1.0 29.4 ± 1.5 25.0 ± 1.5 HET 60meV 21.0 



Tj„/2@ 5K 1.6 ± 0.5 3.0 ± 0.5 4.8 ± 0.8 HET 20/60 meV 

r?„/2@ 5K 1.8 ± 0.5 3.4 ± 0.8 4.7 ± 0.8 HET 60 meV 

rJ„/2@75K 2.5 ± 0.5 3.7 ± 0.5 4.4 ± 0.8 IN6 3.1 meV 

Wo ^(r928 XT97 ^0856 

B40 0.052 0.069 0.063 

1^441 0.128 0.088 0.089 



TABLE II: crystal-field energies Ei,2, widths (HWHM) T/2, and crystal-field parameters are given in meV. The fifth column 
gives instrument and incident energy Ei„ from which the parameters were determined. Values of CeAu2Si2 are shown for 
comparison. The crystal-field parameters from the full multiplet calculation are given in Stevens formalism. 



suiting ground state wave functions, with an arbitrarily 
chosen phase, are: 



CeRhlns: 


|0) 


— ^ 7 


= 0.62| ±5/2) - 


h 0.781 T 3/2) 


Celrlus: 


|0) 


— ^ 7 


= 0.50|±5/2)- 


1-0.871 T 3/2) 


CeCoIus: 


|0) 


- 

— ^ 7 


= 0.36|±5/2)- 


h0.93| t3/2) 



The orbitals shown in Figure [T] show the spatial distri- 
butions of the 4/ electrons for the respective crystal-field 
ground states. The higher the |5/2) contribution to the 



ground state the fiatter the 4/ distribution, i.e. CeRhlns, 
which does not become superconducting at ambient pres- 
sure, has the flattest 4/ orbital. The pure |5/2) orbital is 
donut and the pure |3/2) is yo-yo shaped (see e.g. Willers 
et a].—). The general trend of a decreasing |5/2) con- 
tribution to the ground state from M = Rh, Ir to Co 
agrees with the INS findings by Refs. [29ll30l but we ob- 
serve smaller |5/2) contributions (note: a/(1 — a^) = (3 
as given by Christiansen et al). This is most likely due 
to systematic errors in the phonon correction of the INS 
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neutron data (see appendix). 



B. INS: crystal-field transition energies 

Figure [2] shows the scattering function S{20, cj) at T = 
5 K for smah scattering angles 20— 19° measured with 
two incident energies, Ein = 20 meV a) to e) and 60 meV 
f) to h). The separation of magnetic and phonon correc- 
tion has been performed as described in the appendix. 

= 20 meV: a) and b) exhibit the scattering function 
S{20,uj) of two non-magnetic reference samples and c) 
to d) show S{29, w) of the three cerium compounds. The 
blue triangles are the total scattering. The data of the 
La samples show that there is only little phonon scat- 
tering for 20 meV incident energy. It has nevertheless 
been corrected for (method 2) in order to determine the 
magnetic scattering in the cerium data. The red circles 
are the pure magnetic scattering after phonon correction. 
Ein = 60 mcV: the blue triangles in f) to h) are the total 
scattering for small scattering angles, i.e. for low Q and 
the black lines reflect the phonon scattering as obtained 
from the high to low angle scaling method. The phonons 
scale with R = 1/8 from 29^ 135° to 20= 19° as empiri- 
cally found from the non-magnetic reference samples (see 
appendix). The red circles are the pure magnetic scat- 
tering in the low angle banks, resulting from subtracting 
the scaled phonon intensities (black lines) from the total 
scattering (blue triangles). 

All three compounds exhibit two magnetic ground 
state excitations at about 5-7 and 24-30 meV. The low 
energy one (CFl) is best resolved in the 20 meV data 
while the 60 meV data show both excitations (CFl and 
CF2). The quasielastic scattering cannot be resolved 
with these thermal measurements. The inelastic, mag- 
netic excitations have been described with two inelastic 
Lorentzians. The crystal-field transitions are fairly sharp 
in CeRhIn5 but much broader for Celrlns and CeCoIns, 
leading to considerable error bars in the line widths and 
position. The 60 meV data are described with two inelas- 
tic Lorentzians of same widths, but only the parameters 
of the higher energy excitations were freely varied, posi- 
tion and width of the low energy excitations have been 
determined from the 20 meV data. We can state that the 
excitations broaden from M = Rh to Ir and Co, consis- 
tent with the expected increase in the electronic specific 
heat coefficient 7 w 1/Tk from Rh to Ir to Co. Table 2 
gives the summary of all parameters. 

In the Ein — 60 meV data phonon and magnetic scat- 
tering are very much intermixed, the line widths are 
broad and we know from the non-magnetic reference sam- 
ples that there is a substantial error in not taking out 
enough phonon scattering below 10 meV, i.e. in the range 
of the lower crystal-field excitation (this is explained in 
the appendix, see figure |6]). We know further that the 
fine structure of the phonon scattering is not ideally rep- 
resented by the Q scaled high angle data (see appendix, 
right hand panel of figure [B]). Hence we do not attempt to 



determine magnetic intensities to obtain transition ma- 
trix elements. The magnetic signal is nevertheless strong 
enough to determine the crystal-field energies fairly well, 
irrespective of the method of phonon correction as can be 
seen from the fact that we obtain similar line positions 
as Christianson et al. who applied some extra scaling 
factor for the phonon correction in addition to the one 
found from the La samples (compare appendix with Refs. 
[29II30I ). However, for Celrlns Christianson et al. can only 
state that the low lying crystal-field state must be some 
where between and 7 meV. Our neutron data show, due 
to the better flux and resolution, that the level is located 
between 5.2 ± 1 meV (see figure [2] e)). This is supported 
by the fact that the polarization effect of the XAS data 
increases from 20 to 80 K (see section C, figure [3l), i.e. 
when the thermal occupation of the first excited state 
takes place, before it decreases again due to the begin- 
ning population of the second excited state. 

We find narrower line widths and we find that the 
width of the spectra increases in the sequence M=Rh 
to Ir and becomes broadest for Co in contrast to Refs. 
[29II30I . We believe that this is due to the better signal to 
noise ratio of the present data. While we do not consider 
the quality of the magnetic intensities to be sufficiently 
accurate to determine the wave functions pro perly from 
these data, we agree with Christianson etal. I29ll30l that 
general intensity considerations based on the selection 
rule AJz = ±1 lead to the conclusion that the |5/2) con- 
tribution to the ground state decreases from Rh to Ir and 
Co, that is in agreement with our low temperature XAS 
data. We conclude further from the intensity ratios, as 
Chrstianson etal . [29ll30l . that the second excited state is 
the pure |l/2) in all compounds. 

The same set of data was taken at T = 120 K (not 
shown here) where we would expect to see the transition 
from the first |1) to the second |2) excited state at E2-E1, 
but we do not resolve another peak. However, it should 
be noted that at 120 K the phonon correction is more 
important with respect to 5 K while the magnetic lines 
are broader. In addition, this third transition is weaker 
than the two ground state excitations. 



C. Temperature dependent polarized soft XAS: 
sequence of crystal-field states 

When analyzing the temperature dependence of the 
polarization dependent XAS data we will only show the 
M5 edge for clarity. Figure [3] shows the temperature 
dependence of the LD for all three cerium compounds. 
For each compound the measured data and correspond- 
ing simulations are shown. For CeRhlns the LD in the 
experimental data increases slightly from 18 to 80 K in- 
dicating the population of the first excited crystal-field 
level and then decreases with further rising temperature 
due to population of the next state. When all states are 
equally populated the polarization should vanish entirely 
since a equal occupation resembles the fully degenerate 
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FIG. 3: (color online) Left panel: Temperature-dependent linear polarized XAS spectra at the Ce^+ M5 edge of CeRhlns, 
Celrlns, and CeCoIns. The solid lines correspond to the measured data, the dotted lines to the simulation based on the 
crystal-field energies from the neutron data and the ground state wave functions from the low temperature XAS data. 




Photon energy in eV P'^'""" o"^'^"/ 



FIG. 4: (color online) Isotropic spectra of CeMIus with 
M=Rh, Ir and Co, (at 20 K) and of CeAu2Si2 (at 40 K) 
for comparison. All spectra were recorded under the same 
beam line conditions at the NSRRC. In the right panel are 
blow-up graphs with arrows pointing out the 4/" initial state 
contributions. 



Hund's rule ground state which has spherical symme- 
try. To the right of the measured data the simulated 
XAS spectra are shown. There, the crystal-field energies 
as determined in the neutron experiment and the mix- 
ing factor a as obtained from the low temperature XAS 
data are used as input parameters. The |l/2) has been 
assumed to be the highest lying crystal-field level and 
only the population of states has been adapted to the 
corresponding temperature. The simulation reproduces 
well the trend of the temperature dependence of the LD, 
i.e. it increases at first and then decreases above 80 K 
with further rising temperature. The assumption of a dif- 
ferent order of states would lead to a change of sign in the 
polarization for increasing T, which is in contradiction to 
the observation and can therefore be excluded. The other 
panels of figure[3]show the same set of spectra for Celrlns 
and CeCoIns. For Celrlns we observe very much the 
same effect as for CeRhlns and we are able to simulate 
the general trend of the temperature dependence in the 
same manner. For CeCoIns the temperature dependence 
seems different, but is based on the same effect, namely 
occupation of the next higher states: the low tempera- 
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ture polarization has a different sign with respect to M= 
Rh and Ir (see section A) so that here the occupation of 
the first excited state leads to a decrease of polarization 
with rising temperature and a change of sign between 50 
and 80 K. Then the LD increases at first with further 
rising temperature and starts to decreases again above 
180 K. Here too we are able to simulate this temperature 
trend under the assumption that the [1/2) is the highest 
lying crystal-field state. For M=Co, the assumption of a 
different order of states would not give rise to a change of 
sign in the polarization as function of temperature. This 
analysis of the temperature dependence is analogous to 
our previous results on CePtaSi.'^^ Ground state, energy 
splittings, and order of states describe the crystal-field 
potential fully and the corresponding crystal-field param- 
eters (in Stevens formalism) are listed in table 2. 

While the temperature dependence of the LD is qual- 
itatively reproduced with our simulations, there is also 
some quantitative deviation: at high temperatures the 
measured LD is always smaller than the simulated one. 
We exclude depolarization effects due to a) surface de- 
grading and/or b) polycrystalline contributions. Possibil- 
ity a) can be excluded since we recleaved and repeated the 
measurements to assure reproducibility (see table 1) and 
we exclude b) since it would require an unrealistic 35% 
of polycrystalline contribution in order to account for the 
mismatch in e.g. the Co data at 280 K. More physical 
and interesting is to consider Kondo interactions: they 
will have a depolarizing effect too. The hybridization of 
the 4/^ state with the surrounding conduction band has 
not been considered in the present analysis of the XAS 
data. Yet, the existence of the latter can be seen from 
the 4/° initial state satellites at the high energy tails of 
the M4^5 edges. Figure |4] shows the isotropic spectra 
and the tails of the M4_5 edges on a blown-up scale for 
the CeMIns compounds and in comparison for CeAu2Si2. 
The arrows in Figure |4] indicate the position of the /° 
spectral weight. CeAu2Si2 is an antiferromagnet with a 
small hybridization temperature (T* — 1.5 K)^^ and with 
a good agreement of simulated and measured XAS data 
at all temperatures25.. It is interesting to note that the 
4/° spectral weight is basically non-existent in CeAu2Si2, 
but stronger in the CeRhlns data and again more pro- 
nounced for Celrln5 and CeCoIns. In the same sequence 
the inelastic line widths in the INS spectra and the devi- 
ation between simulated and measured spectra increase. 
We speculate that hybridization effects, which are not yet 
included in the calculation, may be responsible for these 
quantitative discrepancies. It is desirable that further 
theoretical work be carried out, using for instance the 
Anderson impurity model, to find out how much the CF 
wave functions are modifed from our present estimates. 
We nevertheless expect that these modifications are very 
modest for the CeRhIn5 since the /° weight in the ground 
state is minimal. We may even speculate that the correc- 
tions are also small for the CeCoIns and Celrlns in view 
of the special condition that the crystal-field ground state 
wave function is very close to cubic, i.e. almost isotropic. 



D. INS: hybridization temperature T* 

Figure [5] shows high resolution inelastic neutron scat- 
tering data of CeMIns M=Rh, Ir and Co, taken with 
IN6 at ILL with Ei„ = 3.1 meV incident energy at 1.5, 
8, and 75 K. The blue triangles (pointing down) are the 
total scattering from the cerium samples. We further 
show the scattering from Lalrlns at 75 K (see black tri- 
angles pointing up in the 1.5 and 75 K CeRhIn5 spec- 
tra) in order to verify that the phonon scattering is 
negligible at these temperatures in this energy window. 
At 1.5 K CeRhIn5 is magnetically ordered (TAr=3.8 K) 
and the scattering function exhibits in addition to the 
elastic line (incoherent, elastic, nuclear scattering) some 
magnon density of states which peaks at about 1.7 meV. 
At 8 K the spectrum consists mainly of quasielastic 
scattering, which is well described with a quasielastic 
Lorentzian with HWHM r/2 = 1.2 (±0.1)meV. Note, 
that for hT/2 < ksT a quasielastic Lorentzian appears 
highly asymmetric in energy^. The crystal-field excita- 
tion CFl contributes only little in this energy window 
but has nevertheless been taken into account according 
to the crystal-field model (see green line). A quasielastic 
line width of r/2 = 1.2 (±0.1) meV corresponds with hoj 
= ksT to T* « 14 (±1) K. At 1.5 and 8 K the scattering 
intensity for CeCoIns and Celrlns is considerably lower 
with respect to CeRhlns, but not zero, so that we con- 
clude the quasielastic lines are broader and peaked out- 
side the energy window available at low temperatures. 
It should be mentioned that although CeCoIns is in the 
superconducting phase at 1.5 K, we did not observe the 
spin resonance at 0.6 meV— , most likely because this 
resonance is fairly sharp in reciprocal space whereas the 
present data are averaged over all Q. 

Before discussing the CeCoIns and Celrlns data at 1.5 
and 8 K data in more detail, we look at the high res- 
olution spectra at 75 K (right hand panel of figure [S]) . 
Here the energy window is larger due to population of 
states on the neutron energy gain side. Thanks to popu- 
lation and resolution the quasielastic scattering and the 
low lying crystal-field excitation can be observed simulta- 
neously in the same energy window. The black triangles 
in the CeRhlns spectrum are the scattering from Lalrlns 
at 75 K. Since it is almost negligible we will not con- 
sider it when describing the cerium data. For CeRhlns 
quasi- and inelastic scattering are well resolved while for 
Celrlns and CeCoIns the spectra are broader so that it 
is difficult to separate the two. We fit the data with the 
crystal field model, i.e. keep the intensity ratios of quasi- 
and inelastic lines as well as the line positions fixed (see 
parameters in table 2). Thus only the line widths and 
an overall intensity parameter were varied. The latter 
is necessary although we scaled to absolute intensities, 
most likely since the absorption corrections are so large. 
The data are well described with our crystal-field param- 
eters, the resulting quasi- and inelastic line widths at 75 
K are listed in table 2. 

We return to the 8 K data in order to give an estimate 
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FIG. 5: (color online) Neutron data with 3.1 meV incident energy. For T = 1.5 and 8 K the neutron energy loss, for 75 K the 
neutron energy gain side is shown. The blue triangles (pointing down) are the total scattering of the cerium samples, the black 
triangles (pointing up) are the La scattering. The 75 and 8 K data are described with one quasi- (black) and one inelastic 
Lorentzian (green) with intensity ratios according to the crystal- field models (see table 2). The red lines reflect the total, fitted 
magnetic scattering. The line widths are discussed in the text. 



for the quasielastic line widths of CeCoIns and Celrlns. 
The quasielastic intensity at 8 K is scattering from the 
crystal-field ground state which is proportional to the 
square of the transition matrix element. The ground 
state of each compound is know from the low temperature 
XAS data, hence we know the respective matrix elements 
of the quasielastic scattering. We further know the rela- 
tive scaling factor of the data from the 75 K spectra. For 
CeRhlns the quasileastic width and intensity can easily 
be fitted. Knowing scaling factors and matrix elements 
allows us then to estimate the line widths of the Co and 
Ir samples at 8 K (middle columns of figure [5]) . We use 
the widths at 75 K as a crude guess for the quasielastic 
line widths of CeCoIn5 and Celrln5 at 8 K and take into 
account the inelastic scattering according to the crystal- 
field model. The red lines are the result. The agreement 
is reasonable so that we give the following values for the 
corresponding coherence temperatures T*: 



CeRhlns: T* = U±1K 
Celrlng: T* « 30 ± SiiT 
CeCoIns: T* « 45 ± 8i^ 



IV. DISCUSSION 

The hybridization temperatures T* as determined from 
the present neutron scattering experiments are smallest 
for CeRhlus and largest for CeCoIus in agreement with 
the increasing spectral weight of the contribution to 
the XAS spectra. The values we find from neutron scat- 
tering agree well with temperatures below which Fermi 
liquid behaviour sets in: the anomalous Hall effect below 
20 K for CeRhIng, 31 K for Celrlng, 53 K for CeCoInsi". 
Knight shift experiments give 10-12 K for CeRhlns, and 
for CeCoIug 42 K for in plane and 89-95 K for out of 
plane, i.e. a powder averaged value of about 58 K—. Ther- 
mal and transport measurements by Nakatsuji et al. give 
T* 45 K for CeCoIug and the entropy of the specific 
heat of CeRhlns reaches l/2i?ln2 at about 10-12 K4. 
Our findings are further in agreement with the findings 
of the Fermi surfaces. While the general features of the 
Fermi surface of CeRhlns are more like the Fermi sur- 
face of LaRhlns, which has no 4/ electrons, the Fermi 
surfaces of CeCoIus and Celrlus are well described with 
a more itinerant 4/ band model^. A summary of val- 
ues and references can be found in the supplementary 
information of Ref. [111. 
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CeRhlns is the most localized member of this fam- 
ily. The hybridization temperature of CeRhlns com- 
pares with those of other heavy fermion materials like 
CeCu2Si2 (T* « 10 K}^, CeRuaSia (T* « 10 K)^^, 
and with the ones of the antiferromagnetic compounds 
CePd2Si2 (T* = 10 K, = 8 Kp- and the cubic par- 
ent compound Celng (T* w 10 K, T^v = 10 K)^. With the 
exception of CeRuaSia all these compounds, exhibit su- 
perconductivity: CeCuaSia at ambient pressure depend- 
ing on sample stoichiometry, or CeRhlus, CePdaSia and 
Celna with an applied pressure of 1.6 GPa and about 
2.5 GPa for the latter two^i^. It is intriguing that 
for CeAuaSia, which also orders antiferromagnetically at 
Tn ~ 10 K but has a much smaller hybridization tem- 
perature of T* = 1.5 K^*, no superconductivity has been 
reported up to 17 GPa^^. Since pressure on the cerium 
ion increases delocalisation, these findings underline the 
idea that sufhcient Kondo screening favors superconduc- 
tivity to the detriment of antiferromagnetic order^' For 
Celrlus and CeCoIns where the hybridization temper- 
atures are larger the Kondo screening seems to be suf- 
ficiently large so that superconductivity can develop at 
ambient pressure. Here the larger hybridization temper- 
ature of CeCoIus (Tc = 2.3 K) with respect and Celrlns 
(Tc = 0.4 K) goes along with a higher superconducting 
transition temperature. 

There is another aspect which makes Celrlns and 
CeCoIns remarkable. They are rare examples for com- 
pounds where the hybridization temperature T* is sim- 
ilar in magntitude to the size of the (low lying) crystal- 
field splitting. One may speculate that this should have 
an effect on the degeneracies involved when describing 
ground state properties. 

Along with the increasing hybridization from M=Rh, 
via Ir to Co goes a decrease of the |5/2) contribution 
to the ground state. CeRhlns which has the flattest 4/ 
orbital (see orbitals in figure 2) does not become super- 
conducting at ambient pressure and it does not appear in 
the Tc versus c/a scaling plot as suggested by Pagliuso^^, 
although its c/a ratio is in between the values of Celrlns 
and CeCoIns. The CeIri_a;Rha;In5 and CeCoi_3;Rhj;In5 
systems, however, exhibit superconductivity and fit into 
this c/a scaling. It would be interesting to see where, 
in a similar scaling of T^ with the 4/ wave functions, 
the CeIri_a;Rha;In5 and CeCoi_a;Rha;In5 would fit. Since 
the crystal-field energies do not vary much from sample 
to sample it would be sufficient to determine the ground 
state wave functions with linear polarized XAS at low 
temperatures. 



V. SUMMARY 

We have determined the hybridization temperatures 
T* and crystal-field schemes of CeMIns M=Rh, Ir and 
Co with inelastic neutron scattering and polarized soft 
x-ray absorption. The hybridization temperatures T* as 
determined from the line widths of the inelastic neutron 
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FIG. 6: (color online) Inelastic neutron scattering data of 
the non- magnetic reference compounds YCoIns and Lalrlns. 
Left: high and low angle data at 5 K. Right: 5 K low an- 
gle data shown on an expanded intensity scale and scaled 
(7? — 1/8) high angle data. Bottom right: the red squares 
are the same low angle data taken at 280 K. The orange line 
represents the 5 K data scaled to 250 K with the Bose factor. 



data increase from M = Rh to Ir, and are largest for 
Co which supports the idea that increasing Kondo inter- 
action favours superconductivity while preventing long 
range magnetic order. The hybridization temperature 
of CeRhlns, the most localized member of the family, is 
comparable to the CeRuaSia and CeCuaSia. For Celrlns 
and CeCoIns the energy scale of the 4/ conduction elec- 
tron interaction is of the order of the energy of the low 
lying crystal-field excitation, which may have an impact 
on the ground state degeneracy and/or properties. Our 
finding of the crystal-field schemes is coherent with previ- 
ous work by Christianson et al^^^ but we can give more 
precise values for the ground state wave functions from 
our XAS data. We find that the |5/2) contribution to the 
ground state is largest (smallest) for CeRhlns (CeCoIns) 
so that CeRhlns has the flattest 4/ orbital. 



VI. APPENDIX 

A single phonon scattering process increases with mo- 
mentum transfer Q as and magnetic scattering de- 
creases with the magnetic form factor of the magnetic 
ion. Phonon scattering follows Bose and the occupation 
of crystal-field states is according to Boltzmann statis- 
tics. We shall now discuss the two ways of separating 
magnetic and phonon scattering: 1) subtraction of the 
high Q data from the low Q data after scaling the high Q 
data with a scaling factor R which has been determined 
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with a non-magnetic reference compoundiS and 2) sub- 
traction of S{Q,uj) of a non-magnetic reference which 
has been scaled by the average scattering cross-section. 
Method 1): since we group angles and not momentum 
transfers Q we will determine the scaling factor for high 
to low scattering angles 29 and not for Q. The left hand 
panel of figure [5] shows S{29,uj) of YCoIn5 and Lalrlus 
for 60 meV incident energy for large scattering angles 29 
— 135° as black triangles (pointing up) and small scatter- 
ing angles 29 = 19° as blue triangles (pointing down). All 
scattering of these non-magnetic samples is of phonon ori- 
gin so that a phonon scaling factor R from large to small 
scattering angles can be determined. Empirically we find 
for both non-magnetic samples that the high angle in- 
tensities scale to the low angle ones with a factor oi R = 
1/8. The quality of this scaling is shown in the right hand 
panel of figure [6] The blue triangles (pointing down) 
are again the low angle scattering, now shown on an ex- 
panded intensity scale and the black line is the high angle 
scattering scaled by i? = 1/8. While the overall intensity 
is well described, the low angle phonon scattering below 
10 meV is under estimated by such a scaling and the fine 
structure of the low angle phonon scattering above 10 
meV is not so well reproduced. Method 2): next we check 
method 2 by scaling the two non-magnetic reference com- 
pounds to each other with the averaged nuclear cross- 
section, again for grouped angles. The ratio of the aver- 
aged nuclear cross-sections of YCoIns and Lalrlns is 1.37. 
For large scattering angles the 60 meV data seem to scale 
rather well with this value; we find a scaling factor of 1.3. 
However, for small scattering angles they do not. Here 
we find that the two non-magnetic data sets scale best 
with a factor of 1.1. This deviation is probably due to 
multiple scattering which seems stronger in the forward 
detectors. This finding makes this type of scaling some 
what arbitrary for the magnetic samples unless a detailed 



phonon simulation and multiple scattering calculation is 
performed. Another obstacle of method 2) can be that 
often the non-magnetic reference samples are only mea- 
sured at base temperature and the higher temperatures 
are obtained from Bose scaling. The bottom right panel 
of figure [Ushows that this can be traitorous: YCoIus was 
measured at 5 (triangles) and 250 K (circles). Scaling the 
scattering function S{29,ll!,T) at 5 K to 250 K with the 
Bose factor [S{9,uj,250K) = B{250K)/B{5)S{Q,u},5K) 
and n{u},T) + 1 = B{T) = 1/(1 - exp{-huj /keT)) ] 
gives the orange line which highly overestimates the ac- 
tual scattering at 250 K. The Bose scaling does of course 
not take into account multiple phonon processes. 

We therefore discard the direct subtraction method for 
the phonon correction of the 60 meV data and rather ap- 
ply the high to low Q or large to small scattering angle 29 
scaling. However, for the 20 meV data which exhibit very 
little phonon scattering in the forward detectors (see left 
column of figure [2|) the cross-section scaling looks fine 
and we correct the 20 meV data by subtracting the data 
of the cross-section scaled non-magnetic reference sam- 
ples since the high angle data of the 20 meV data contain 
still a non-negligible amount of magnetic scattering due 
to the smaller Q values. 
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